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Abstract 



The shadowing and antishadowing of nuclear structure functions in the Gribov- 
Glauber picture is due respectively to the destructive and constructive interference 
of amplitudes arising from the multiple-scattering of quarks in the nucleus. The 
effective quark-nucleon scattering amplitude includes Pomeron and Odderon contri- 
butions from multi-gluon exchange as well as Reggeon quark-exchange contributions. 
We show that the coherence of these multiscattering nuclear processes leads to shad- 
owing and antishadowing of the electromagnetic nuclear structure functions in agree- 
ment with measurements. This picture leads to substantially different antishadowing 
for charged and neutral current reactions, thus affecting the extraction of the weak- 
mixing angle Ow- We find that part of the anomalous NuTeV result for 9 W could be 
due to the nonuniversality of nuclear antishadowing for charged and neutral currents. 
Detailed measurements of the nuclear dependence of individual quark structure func- 
tions are thus needed to establish the distinctive phenomenology of shadowing and 
antishadowing and to make the NuTeV results definitive. 
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1 Introduction 



The precise determination of the weak-mixing angle sin 2 9w plays a crucial role in 
testing the standard model of electroweak interactions. Until recently, a consistent 
value was obtained from all of the electroweak observables pQ. However, the NuTeV 
Collaboration [2] has determined a value for sin 2 6\y from measurements of the ratio 
of charged and neutral current deep inelastic neutrino-nucleus and anti-neutrino- 
nucleus scattering in iron targets which has a 3a deviation with respect to the fit 
of the standard model predictions from other electroweak measurements pQ. This 
contrasts with the recent determination of sin 2 6\y from parity violation in Moller 
scattering which is consistent with the standard model [2]. Although the NuTeV 
analysis takes into account many sources of systematic errors, there still remains the 
question of whether the reported deviation could be accounted for by QCD effects 
such as the asymmetry of the strange-antistrange quark sea pfl |H] or other Standard 
Model effects |H1 CZl [HI El Hill UTJ [T21 UHl [131 EISl [TE] - In this paper we shall investigate 
whether the anomalous NuTeV result for sin 2 6\y could be due to the different behavior 
of leading-twist nuclear shadowing and antishadowing effects for charged and neutral 
currents. 

The physics of the nuclear shadowing in deep inelastic scattering can be most 
easily understood in the laboratory frame using the Glauber-Gribov picture [T7| I18j. 
The virtual photon, W or Z° produces a quark- ant iquark color-dipole pair which can 
interact diffractively or inelastically on the nucleons in the nucleus. The destruc- 
tive interference of diffractive amplitudes from pomeron exchange on the upstream 
nucleons then causes shadowing of the virtual photon interactions on the back-face 
nucleons [THl EI3 ED E21 EH1 E3j • As emphasized by Ioffe |21j, the coherence between 
processes which occur on different nucleons at separation La requires small Bjorken 
xb '■ 1/Mxb = 2v/Q 2 > La- The coherence between different quark processes is also 
the basis of saturation phenomena in DIS and other hard QCD reactions at small 
xb [25] . and coherent multiple parton scattering has been used in the analysis of 
p + A collisions in terms of the perturbative QCD factorization approach 26J. An 
example of the interference of one- and two-step processes in deep inelastic lepton- 
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nucleus scattering illustrated in Fig. [T] 




Figure 1: The one-step and two-step processes in DIS on a nucleus. If the scattering 
on nucleon Ni is via pomeron exchange, the one-step and two-step amplitudes are 
opposite in phase, thus diminishing the q flux reaching N%. This causes shadowing of 
the charged and neutral current nuclear structure functions. 

An important aspect of the shadowing phenomenon is that the diffractive con- 
tribution 7*iV — > XN' to deep inelastic scattering (DDIS) where the nucleon N\ in 
Fig. ^ remains intact is a constant fraction of the total DIS rate, confirming that it 
is a leading-twist contribution. The Bjorken scaling of DDIS has been observed at 
HERA [23EH1E2I- As shown in Ref. |U|, the leading- twist contribution to DDIS 
arises in QCD in the usual parton model frame when one includes the nearly instan- 
taneous gluon exchange final-state interactions of the struck quark with the target 
spectators. The same final state interactions also lead to leading-twist single-spin 
asymmetries in semi-inclusive DIS [3T]. Thus the shadowing of nuclear structure 
functions is also a leading-twist effect. 
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It was shown in Ref. [32] that if one allows for Reggeon exchanges which leave 
a nucleon intact, then one can obtain constructive interference among the multi- 
scattering amplitudes in the nucleus. A Bjorken-scaling contribution to DDIS from 
Reggeon exchange has in fact also been observed at HERA [2*%| I2TJ] . The strength 
and energy dependence of the C = + Reggeon t— channel exchange contributions 
to virtual Compton scattering is constrained by the Kuti-Weisskopf [33] behavior 
F2(x) ~ of the non-singlet electromagnetic structure functions at small x. 

The phase of the Reggeon exchange amplitude is determined by its signature factor. 
Because of this complex phase structure [32J, one obtains constructive interference 
and antishadowing of the nuclear structure functions in the range 0.1 < x < 0.2 - a 
pronounced excess of the nuclear cross section with respect to nucleon additivity |33j . 

In the case where the diffractive amplitude on N\ is imaginary, the two-step 
process has the phase i x i = —1 relative to the one-step amplitude, producing 
destructive interference. (The second factor of i arises from integration over the 
quasi-real intermediate state.) In the case where the diffractive amplitude on Ni 
is due to C = + Reggeon exchange with intercept or(0) = 1/2, for example, the 
phase of the two-step amplitude is ^ (1 — i) x % = -^(i + 1) relative to the one- 
step amplitude, thus producing constructive interference and antishadowing. This is 
discussed in more detail in the following sections. 

Odderon exchange due to three-gluon exchange leads to an elastic quark-nucleon 
amplitude which is nearly real in phase, thus providing an additional mechanism 
for antishadowing. We shall show that the combination of Pomeron, Reggeon, and 
Odderon exchanges in multi-step processes leads to shadowing and antishadowing of 
the electromagnetic nuclear structure functions in agreement with measurement in 
electromagnetic interactions. Shadowing of the nuclear structure functions is thus 
due to the dynamics of 7* A interactions; it is not a property of the nuclear light-front 
wavefunction computed in isolation [30] . 

Evidence for the Odderon has been illusive; a detailed discussion can be found in a 
recent review by Ewerz [35J. A clear signal appears in the difference of proton-proton 
vs. proton anti-proton scattering. From the perspective of QCD, the Odderon repre- 
sents the color-singlet effects of three gluons in the t— channel. A general treatment in 
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the context of the BFKL program has been given by Bartels, Lipatov, and Vacca [5o] 
The Odderon has Regge intercept «o ~ 1, C — — , and thus its phase is nearly pure 
real. The Odderon does not contribute directly to the structure functions since it 
gives a real contribution to the virtual Compton amplitude. However, it can play an 
important role in the multi-scattering series in the nuclear target. 

There can be other important antishadowing mechanisms. Processes which can 
occur on a nucleus, but are forbidden on a nucleon, will enhance the nuclear structure 
functions. For example, pseudoscalar Reggeon exchange amplitudes do not contribute 
to DIS on a nucleon target since the helicity-conserving forward amplitude 7*iV — > 
7*iV vanishes at t — 0. However in the nuclear case, the interactions of the scattered 
quark (due to pomeron exchange) on a second nucleon N 2 in a nuclear target can skew 
the kinematics of 7*iVi — > 7*iV{, thus allowing the pseudoscalar exchange to occur on 
the nucleon Ni at t 7^ 0. This also requires nonzero orbital angular momentum of the 
nucleons in the nuclear wavefunction. Notice that the virtual Compton amplitude on 
the nucleus 7* A — > 7* A is still evaluated at zero momentum transfer t = 0. Thus in 
general one should include pseudoscalar exchange in the parametrization of the quark 
multiple scattering processes. 

In Fig.|2]we illustrate several leading- twist QCD contributions to the nuclear struc- 
ture function as calculated from the absorptive part of the forward virtual Compton 
amplitude Im T{^*A —>■ 7* A), in the q + — 0, q\ — Q 2 parton model frame and in the 
laboratory frame where q + > 0. Notice the final-state two-gluon exchange "pomeron" 
interaction of the outgoing quark on a target neutron. Figure 121 (c) is an illustration of 
a doubly inelastic discontinuity of the same two-step process as (a) in the laboratory 
frame q + > 0. The uu fluctuation of the virtual photon first scatters inelastically on 
a neutron via a single gluon exchange which produces an excited color state 8c of 
the neutron. This is then followed by the annihilation of the u quark on a proton. 
The two-step amplitudes of (b) or (c) will interfere destructively with the single-step 
annihilation amplitude on the proton alone, thus producing shadowing. If the pro- 
ton spin S% is flipped (A' 7^ A) by the valence interaction, as occurs in pseudoscalar 
Reggeon exchange, then the single-step process cannot contribute to the forward vir- 
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Figure 2: Representation of leading-twist QCD contributions to the nuclear struc- 
ture function as calculated from the absorptive part of the forward virtual Compton 
amplitude Im T(7*v4 — > j*A). (a) Illustration of a two-step contribution in the usual 
q + = 0, q\ = Q 2 parton model frame. The deep inelastic scattering of a lepton on 
a valence quark of a target proton is followed by the final-state two-gluon exchange 
"pomeron" interaction of the outgoing quark on a neutron, (b) Illustration of the 
physics of the same two-step process illustrated in (a), but in the laboratory frame 
where q + > 0. The uu fluctuation of the virtual photon scatters elastically via two- 
gluon exchange on a neutron; this is then followed by the annihilation of the u quark 
on a proton, (c) Illustration of a doubly inelastic discontinuity of the same two-step 
process as (a) in the laboratory frame q + > 0. The uu fluctuation of the virtual pho- 
ton first scatters inelastically on a neutron via a single-gluon exchange which produces 
an excited color state 8c of the neutron. This is then followed by the annihilation of 
the u quark on a proton. 
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tual Compton amplitude, and the two-step process itself produces antishadowing of 
the valence quark distributions. Similar processes occur in the case of the electroweak 
currents. 

Figure El illustrates a similar situation, but for the three-gluon "Odderon" ex- 
change. In this case, the two-step amplitudes of (b) or (c) can interfere constructively 
with the Regge-behaved single-step annihilation amplitude on the proton alone, thus 
producing antishadowing. Similar processes occur in the case of the electroweak cur- 
rents. 

The Reggeon contributions to the quark scattering amplitudes depend specifically 
on the quark flavor; for example the isovector Regge trajectories couple differently to u 
and d quarks. The s and s couple to yet different Reggeons. This implies distinct anti- 
shadowing effects for each quark and antiquark component of the nuclear structure 
function; this in turn implies nonuniversality of antishadowing of the charged, neutral, 
and electromagnetic current. Anti- neutrino and neutrino reactions will also have 
different antishadowing effects. In addition, there is another source of antishadowing, 
specific to non-abelian theories, which is discussed in more detail in the appendix. It 
includes one-gluon exchange x Reggeon exchange, assuming the existence of hidden- 
color components in the nuclear wavefunction. We do not explicitly include this effect 
in our analysis since the parameterization is uncertain. 

There are also antishadowing contributions arising from two-step processes involv- 
ing Reggeon x Reggeon exchange, but these contributions are power-law suppressed 
in the Bjorken limit. We will not include these higher-twist effects in our analysis. 

In this paper we shall show in detail that the Gribov-Glauber picture for nuclear 
deep inelastic scattering leads to substantially different nuclear effects for charged 
and neutral currents; in fact, the neutrino and antineutrino cross sections are each 
modified in substantially different ways due to the various allowed Regge exchanges. 
This non-universality of nuclear effects will modify the extraction of the weak-mixing 
angle sin 2 9yy obtained from the ratio of charged and neutral current deep inelastic 
neutrino-nucleus scattering. 
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Figure 3: Representation of leading-twist QCD contributions to the nuclear struc- 
ture function from the absorptive part of the forward virtual Compton amplitude 
Im T(^y*A — » 'j* A) involving odderon exchange, (a) Illustration of a two-step contri- 
bution in the q + = 0, q\ = Q 2 parton model frame - deep inelastic lepton scattering 
on a valence quark of a target proton followed by the final-state three-gluon exchange 
"Odderon" interaction of the outgoing quark on a target neutron, (b) Illustration of 
the physics of the same two-step process shown in (a), but in the laboratory frame 
where q + > 0. The uu fluctuation of the virtual photon first scatters elastically via 
three-gluon exchange on a neutron; this is then followed by the annihilation of the 
u quark on a proton, (c) Illustration of a doubly inelastic discontinuity of the same 
two-step process as (a) in the laboratory frame q + > 0. The uu fluctuation of the 
virtual photon first scatters inelastically on a neutron via two-gluon exchange which 
produces an excited color state 8c of the neutron. This is then followed by the an- 
nihilation of the u quark on a proton. The two-step amplitudes of (b) or (c) can 
interfere constructively with the Regge-behaved single-step annihilation amplitude 
on the proton alone, thus producing antishadowing. Similar processes occur in the 
case of the electroweak currents. 
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2 Nuclear shadowing and antishadowing effects due 
to multiple scattering 



In this section, we will extend the analysis of Ref. jH2| for the electromagnetic interac- 
tion case to the neutrino DIS case. The general approach is based on the "covariant 
parton model [37J 13%] - which provides a relationship of deep inelastic cross section 
to quark-nucleon scattering. The central idea is the following: In neutrino DIS on 
a nucleus A, although the virtual current may interact inelastically with a nucleon 
coming from the nucleus in a one-step process as shown in Fig. ^a), it can also in- 
teract elastically with several nucleons before the final nucleon interacts inelastically 
as depicted in Fig. E^b). The interacting antiquark or quark is spacelike: 

r 2 = -k 2 = x(s + k 2 ± )/(l - x) - xM 2 + k\ (1) 

is the negative of the invariant momentum squared of the interacting parton. Here 
M = \{M p + M n ) is the nucleon mass, and k± is the parton's transverse momentum. 
The quark-nucleon amplitude is assumed to be damped at large quark virtuality r, 
so that the quark-nucleon invariant s = (k + p) 2 grows as 1/x. This description of 
deep inelastic scattering is consistent with recent analysis of the role of final state 
gluon interactions in QCD when one chooses light-cone gauge to make the Wilson 
line integral vanish [3D] . 

At high energies the phase of the elastic amplitude is approximately imaginary 
since it corresponds to Pomeron exchange. The accumulated phase in multiple scat- 
tering is also imaginary. Therefore, the two-step amplitude is coherent and opposite 
in phase to the one-step amplitude where the beam interacts directly on N 2 without 
initial-state interactions; the target nucleon N 2 feels less incoming flux, which results 
in nuclear shadowing. Since there may be an Reggeon or Odderon contribution to 
the qN amplitude, the real phase introduced by such contributions leads to antishad- 
owing effect. In this picture, antishadowing is attributed to a dynamical mechanism 
rather than being enforced to satisfy the momentum sum rule jSH] • The sum rule can 
still be maintained by the nuclear modifications of the gluon distribution. 

We now develop the detailed formulas which describe the nuclear shadowing and 
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antishadowing effects. 



2.1 Parameterizations of quark-nucleon scattering 



We shall assume that the high-energy antiquark-nucleon scattering amplitude Tq^ 
has the Regge and analytic behavior characteristic of normal hadronic amplitudes. 
Following the model of Ref. |S2], we consider a standard Reggeon at oir = |, an 
Odderon exchange term, a pseudoscalar exchange term, and a term at = — 1, in 
addition to the Pomeron-exchange term. 

The Pomeron exchange has the intercept ap = 1 + 5. For the amputated q — N 
amplitude T^v and q — N amplitude T q n with q = u, and d, N = p, and n, we 
assume the following parameterizations, including terms which represent pseudoscalar 
Reggeon exchange. Then resulting amplitudes are: 



T- 

± ti- 



er 



s I / + tan — J /3i(t 2 ) - s(3 (t 2 ) - (1 - i)s 1 ' 2 ^ 2 {r 2 ) 



(2) 
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T 
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s ^ + tan /9l(r 2 ) + s/3o(r 2) _ (1 _ i)a l/2^+ (T 2) (g) 

(1 + z)s 1/2 A ; 2 (^) + +^(1 - i)s 1/2 /3f^ do (r 2 ) 



Hl-^ 2 pi; 2 (T 2 ) + (l + ^ 2 f3lJ 2 (T 2 ) 



1/2 V 



rf-n — c 



8 I i + tan y j A(r 2 ) + S p {r 2 ) - (1 - i^fi^f) 



1/2 V ; (9) 
-(1 + ^^yr 2 ) - (1 - O^^Cr 2 ) + W(l - z) S 1/2 ^/T d °(r 2 ) 

-(l+Oa^Cr 2 ) 

= and 1, since the pseudoscalar term cannot act just once in the multiple 
scattering. 
Here 

-2\ fj 



(10) 



1 + {T 2 /u 2 p 

with j = 1, 1/2, —1, O and pseudo (pseudoscalar). The parameters rij (j = 1, 1/2 
and —1) are taken to be the same as those in Ref. [32] • [See also Table I.] 
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The odd-C Odderon with ao — 1 has a real coupling compared to the imaginary 
coupling of the even C Pomeron. It reduces nuclear shadowing and produces anti- 
shadowing although it does not contribute to the free nucleon structure functions. 
In the following numerical estimate, we take fo — 0.1. In order to fit the large x 
experimental data on the parton distributions of the nucleon, we introduce different 
values for the parameters V 2 _ x which control the off-shell dependence of the q — N 
amplitudes. We denote them as V^J 2 and 77^ 2 for the u and d quarks, respectively. 

We take the overall amplitude strength a to be the same in all cases, with a value: 

u = 66 mb. (11) 

The 1 = 1 Reggeon terms in the amplitudes play a very important role, reflecting the 
sea asymmetry d — u of the nucleon in the low x region. 

In principle, 1 = 1 pseudoscalar exchange should also contribute here, but the 
7 * A — > 7 * A cross section is not sensitive to its parameters, and therefore its 
strength cannot be fixed. Careful fits to DVCS and other processes sensitive to the 
1=1 pseudoscalar exchanges are needed. The 1 = coupling is constrained by our 
fit to antishadowing for electromagnetic DIS. Then we can predict antishadowing for 
weak DIS. 

In principle, each Reggeon in the model of qN scattering amplitude should couple 
to the individual quarks with the appropriate isospin and charge conjugation depen- 
dence. For example, the p Reggeon couples as an I — 1, C — — exchange. Although 
Reggeons of both C — ± appear in the quark-nucleon amplitude, in the end after 
multiple interactions and summing over quark and antiquark currents, the nuclear 
Compton amplitude has only C = + exchange in the t-channel. In our model, the 
Reggeon term in the q — N amplitude with 1 = and C = + is taken to represent 
the sum of the possible Reggeon exchanges. At leading twist only the Pomeron and 
Odderon which derive from gluon exchange survive in the multiple scattering. The 
Reggeon exchange to elastic scattering in the multiple scattering is suppressed. 

In the present analysis, we also include the strange quark contribution. The anti- 
strange quark can scatter elastically on one nucleon via Pomeron, Odderon and <fi 
Reggeon exchanges. The Reggeon intercept for the trajectory is close to a R (0) ~ 
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since a^(m|) = 1 and the Regge slope is universal. Actually, the (ft trajectory can be 
parameterized as an(t) = 0.1 + 0.9t [10]. Then we can parameterize the amputated 
s — N and s — N amplitudes as: 

Tg-N = <Ts-n[™(3{ s) (t 2 ) - s(3^(r 2 ) + s 0A ((l + cosO.Itt) - i sin 0.1tt)Au(t 2 )], (12) 

T S ^ N = a s ^ N [is{3{ s \r 2 ) + s(3 { ^(t 2 ) + s 0A ((l + cosO.lvr) - i sin 0.1tt)(3 .i(t 2 )}, (13) 

where <j-s-n = &s-n = <J and N = p, and n. Since the Pomeron coupling to the 
strange quark could be less in strength than its coupling to light quarks, should 
be smaller than fx — 1 for u and d quarks. We take f[ s ^ = 0.1. The value of the 
Reggeon coupling / .i in Po.iij 2 ) is taken as 0.2 GeV a9 , which is the suitable mass 
dimension in order to have the proper mass dimension of the Reggeon terms. With the 
above choice of the parameters, we can produce a shape of strange quark distributions 
which is close to those obtained by a fit analysis (Hj, CTEQ-5 parametrization ^2] 
and MRST parametrization 43J. 

If we use N = p and n, indicating a proton and neutron target, respectively, then 
for an isoscalar target N , the amplitude T^ No per nucleon is 

T- qNo = l -(T- qp + T- qn ). (14) 



Then we introduce 



and 



T Na (s,T>) = V s (s,r 2 )Aj.(T 2 ), (15) 



iA, (r 2 ) ~ (16) 

Now let us turn to the scattering on a nuclear (A) target. We expect that the 
q — A scattering amplitude can be obtained from the q — N amplitude according to 
Glauber's theory as follows, 

z n i ( Z + N \ '\ 

2W=EE T—TT , , m^ 2 ' 1 ^ (T - )fc2 ^ + k2 ~ l) (17) 
k 1= o fc 2 =o fc i + fc 2 I k\ + k 2 ) M 



where 



M = Min{ki +k 2 ,Z} -Max{ki +k 2 -N,0} + 1 
= Min{kx +k 2 ,N} - Max{kx +k 2 - Z,0} + 1 
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(19) 



47Tp c . m y/2(#2 + 26) 



with 



Pcm. = yjr 2 + (s - M 2 - r 2 ) 2 /4 S , 

i? 2 = -i?o,^o = 1.123A 1/3 fm, 
3 



(20) 



(21) 



and b = 10 (GeV/c) 2 is used [22] • Furthermore, we introduce 



T A (s,r 2 ) = T- qA (s,r 




(22) 



Similar expressions hold for T qA . 

The Regge contribution to the deep inelastic cross section comes from the hand- 
bag contribution to the forward virtual Compton amplitude — > 7*p. The Regge 
behavior x~ aR ^ arises from the summation over higher Fock states. The phase of 
the 1 = Reggeon contribution (— i + with = 1/2 entering the virtual Compton 
amplitude is opposite to the positive imaginary contribution of pomeron exchange 
and thus tends to reduce the deep inelastic cross section on a nucleon. 

As shown in Ref. j3U], the multiple scattering contributions from elastic scattering 
from Reggeon exchange is a higher-twist contribution to the deep inelastic cross sec- 
tion; only gauge interactions have a FSI effect in the Bjorken limit jSHj- The Pomeron 
and Odderon nominally have ao ~ 1, so their contributions to elastic scattering are 
not suppressed in the Bjorken limit, since they are derived from multiple gluon ex- 
change. Thus the Pomeron and Odderon can act any number of times in the nucleus, 
but the Reggeon can act only once at leading twist. In effect the Reggeon does not 
have enough time to form in the FSI at small invariant separation x 2 ~ 1/Q 2 - Thus 
FSIs from Reggeons in T qA with ~ 0.5 should be suppressed in the Bjorken limit 
by a power of 1/Q. In order to implement this we put a suppression factor Ftp in the 
multiple scattering Reggeon terms: 




(23) 



for the u and d quarks. And 




(24) 
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for the s quark with Ql ~ 1 GeV 2 , a typical hadronic scale. 

When we take the limit of large Q 2 , the antishadowing due to elastic Reggeon 
exchanges is quenched; however, the presence of the Odderon can produce anti- 
shadowing. For example, a two-step nuclear process shown in Fig. 3(b) from elastic 
Odderon scattering plus inelastic Reggeon scattering gives a contribution to the vir- 
tual Compton amplitude 1 x i x (i + 1) = (—1 + i). [The middle factor of i is due 
to the cut between the two steps.] The positive imaginary contribution to the two 
step amplitude produces an enhancement of the nuclear cross section relative to the 
nucleon cross section in the regime x ~ 0.1 where the Reggeon contribution to deep 
inelastic scattering is important. This is a key feature of our model. Note that the 
two-step process of Odderon plus Pomeron produces only a real contribution to the 
virtual Compton amplitude. The Pomeron-Pomeron and Pomeron-Reggeon two-step 
contributions reduce the one-step Pomeron plus Reggeon contributions, respectively, 
and thus only produce nuclear shadowing. 

We can also consider the two-step [PiR n ] contribution to Im T y *A- n *A, which 
involves the imaginary part Im R n of the (nonforward) pion Reggeon exchange am- 
plitude. The pion pole term alone is real so we consider the pion Regge trajectory — the 
Reggeized version of pion exchange. When we take its absorptive part, we look at 
the cut through a qq ladder exchanged on the second nucleon. Since the Reggeized 
pion exchange has 1 = 1 and the pomeron is / = 0, the two-step [PiR n ] contribution 
does not contribute to Im T^a-^^a if the target is / = Thus there is no antishad- 
owing contribution from [PiR n ] to the deuteron structure function. However, there 
are other pseudoscalar exchanges possible, such as the 77 Reggeon. 

It is interesting to analyze the situation form the point of view of angular mo- 
mentum. The question is whether the two step process [PiR n ] require orbital angular 
momentum in the ground state nuclear wavefunction. Consider the [PiR n ] contribu- 
tion to the forward virtual Compton amplitude Im T 7 *a-> 7 *a- The pomeron exchange 
on the first nucleon gives a transverse momentum kick kj_. The pion Reggeon ex- 
change on the second nucleon gives a balancing opposite kick — k± so that we can 
have a forward nuclear amplitude. The pseudoscalar exchange on the second nucleon 
is a AL = 1 transition. That is why the amplitude requires nonzero k±. Thus we 
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are actually looking at the overlap of nuclear Fock components of the nucleus with 
AL = 1. This is the same admixture which in the spin-1/2 case which gives a nuclear 
magnetic anomalous moment. 

The unpolarized quark distribution functions in an isoscalar target (No) and nu- 
cleus target (A) are, respectively, 



xq No (x) 



xq (x) 



2 Cx 2 r 

— — / dsd 2 k ± ImT No (s,/i 2 ), 

(2ny i 1 — x J 

(2l)n- x J dsd2k±ImTA{s ^ 2) > 

-T 2 . 



(25) 

(26) 
(27) 



The constant C is related to the parton wave function renormalization constant. 

With the obtained quark distributions xq No( - A \x) for an isoscalar target No or 
a nucleus target A, we can calculate the structure functions for various current ex- 
changes. 

(1) The photon exchange case 



F 1 



7 N (A) 



-\- 

2| 9 



u ( x fo(A) +T[ N (A) 



+ 



d(x) 



N (A) + 2 No{A) 



(28) 



F 7iVo(A) = 2xF lN 0{ A)_ 



(2) The neutral current exchange case 

The structure functions of the NC reaction are 

2 



(29) 



ZN (A) 



(9vf + (9 U A?\ (u N ^ A \x)+u N ^ A \x)) 



+ [(9vf + (9a) 2 ] id NoiA) (x) +d No(A) (x) + s N o(A) ix)+ ^N (A) 



X 



F ZN (A) = 2xF ZNo(A)^ 



(30) 



(31) 



^ZN (A) 



= 2 



g u v g U A(u No{A) (x)-u N ^ A \x)) 



+9v9A(d Na(A) (x) - d No{A) (x) + s N °W(x) - s N °^(x)) 



(32) 
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In the SM the vector and axial-vector quark couplings are given by 

9v - 2 _ 3 sm Uw > 9v ~ ~2 3 A 2' _ 2' 

where sin 2 is the weak-mixing angle. 
(3) The charged current exchange case 
The structure function of the CC reaction is given by 

fY +No{a) = ^^^(IK.P + IK.H+^^^IK.I^K-x) (33) 
+d N ^ A \x)\V ud \ 2 + d N ^ A \Q\V cd \ 2 9(x c - x) 
+s N ^ A \x)\V us \ 2 + s N ^ A \Q\V cs \ 2 6(x c - x), 
here V^- are the Cabibbo-Kobayashi-Maskawa mixing matrix elements. The variable 

x(i+§y (k=c,b), 

x, {k = u, d, s), 



6 



A' 



and the step functions 8(x c — x), 8(xb — x) take into account rescaling due to heavy 
quark production thresholds. 

The structure functions N °( A ) anc i N o( A ) are obtained from (}3*3*j) by the 
replacement of the quark distribution functions q(x, Q 2 ) indicated in the curly brack- 
ets: 

F™ +No(A \x,Q 2 ) = Fr NoiA \x,Q 2 ){q No{A \x,Q 2 )^2x q N ^ A \x,Q 2 ), (34) 

q No{A H^Q 2 )^2Z k q No{A H^Q 2 )}, 
f w+n o( a)^ q2) = 2 F^ +No{A \x,Q 2 ){q N ^ A \x,Q 2 ) ^ -q N ^ A \x,Q 2 )}. (35) 

There are similar formulas for the ^"-current exchange reaction. 
2.2 The values of the parameters 

In the last section we presented the formulas involved in our formalism. We shall take 
the value for most of the parameters to be the same as those in Ref. [32J. The values 
of other parameters are chosen in order to fit the experimental data O HH HE] on 



F 2 °,(F 2 — Fg) , F 2 /F 2 ; they are then checked against the known nuclear shadowing 
and antishadowing effects jlTJ HH1- A summary of the set of parameters is given in 
Tabled 



Table 1: Parameters in our numerical calculation 
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With the above parameters, the average nucleon structure function 

F 2 = = ?L±1L (36) 

for the photon exchange case are shown as a solid curve in Fig. 0] The valence and 
sea contributions to F 2 are also presented as dashed and dotted curves in Fig. F 2 is 
close to the SLAC and NMC jlHj experimental data. We also show our results of the 
difference F 2 — F 2 and the ratio F 2 /F 2 of the nucleon structure functions in Fig. El 
and Fig. |Hl respectively. In the coming subsection we will show the nuclear effects on 
the structure functions and in the next section estimate the nuclear shadowing and 
antishadowing effects on the extraction of sin 2 6 W . 

2.3 Nuclear shadowing and antishadowing effects 

We introduce the ratio 

R = Ff/FJ? (37) 
19 




Figure 4: The calculated nucleon structure function F 2 (biggest solid curve), and 
valence (next solid curve at large x values) and sea (smallest solid curve at large x 
values) contributions to F2, and the corresponding result of the MRST parametriza- 
tion jlH] (dashed curves), at Q 2 = 1 GeV 2 . The experimental data are taken from 
Ref. " 



to indicate the nuclear electromagnetic shadowing and antishadowing effect. We will 
focus on the nucleus 56 Fe since an iron target was used in the NuTeV experiment and 
test the nuclear effect in the x > 0.01 region since 97% of the NuTeV data is from 
0.01 < x < 0.75 [50J. In Figs.[?HEl we show the quark q and anti-quark q contributions 
to the ratio of the structure functions. In order to stress the individual contribution 
of quarks, the numerator of the ratio F A / F^° shown in these two figures is obtained 
from the denominator by a replacement q N °(q N °) into q A (q A ) for only the considered 
quark (anti-quark). Because the strange quark distribution is much smaller than u 
and d quark distributions, the strange quark contribution to the ratio is very close to 
1 although s A /s N ° may significantly deviate from 1. 

In Fig. El we give our prediction of the nuclear shadowing and antishadowing 
effects (the sum of all quark and antiquark contributions to the ratio R) for nuclei 
56 Fe and 40 Ca. From Fig. El we find that our model can explain well the experimental 
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X 



Figure 5: The calculated difference F% — F£ of the nucleon structure functions (solid 
curve), and the corresponding result of the MRST parametrization (dashed curve), 
at Q 2 = 1 GeV 2 . The experimental data are taken from Ref. 




Figure 6: The calculated ratio F^jF^ of the nucleon structure functions (solid curve), 
and the corresponding result of the MRST parametrization [13] at energy scale 1 GeV 
(dashed curve). The experimental data are taken from Ref. 
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Figure 7: The quark contributions to the ratios of structure functions at Q 2 = 1 GeV 2 . 
The solid, dashed and dotted curves correspond to the u, d and s quark contributions, 
respectively. This corresponds in our model to the nuclear dependence of the a(u— A), 
a(d — A), <j(s — A) cross sections, respectively. In order to stress the individual 
contribution of quarks, the numerator of the ratio FA/F^° shown in these two figures 
is obtained from the denominator by a replacement g^ into q A for only the considered 
quark. As a result, the effect of antishadowing appears diminished. 

data on the nuclear shadowing and antishadowing effect in DIS for electromagnetic 
currents. 

We can further check the our model by predicting the ratio of F 2 structure func- 
tions F2A [neutrino] /( 18/5) F2A[muon], which has been measured by the NuTev col- 
laboration The results are shown in Fig lTUl and they agree very well with the 
experimental data, in a calculation with no further free parameters. Notice that the 
data for the ratio tends to go below 1 for x > 0.4, which is also what we find. In 
this case we have taken Q 2 = 20 GeV 2 , which is the average value of the NuTeV 
experiment. 
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Figure 8: The anti-quark contributions to ratios of the structure_ functions at 
Q 2 = 1 GeV 2 . The solid, dashed and dotted curves correspond to u, d and s quark 
contributions, respectively. This corresponds in our model to the nuclear dependence 
of the a{u — A), <j(d — A), cr(s — A) cross sections, respectively. In order to stress the 
individual contribution of quarks, the numerator of the ratio / shown in these 
two figures is obtained from the denominator by a replacement q N ° into q for only 
the considered anti-quark. 



We emphasize that the nuclear shadowing and antishadowing of the different 
currents are not universal since they depend on the different quark species. We still 
have factorization in the sense that we will have the same shadowing quark by quark 
in nuclear Drell-Yan processes. 

In the case of weak currents, the shadowing/ antishadowing effects are strongly 
influenced by the behavior of the structure function F3, which is not present in 
the electromagnetic case. We will present in the next section cross section ratios 
(nucleus/nucleon) to illustrate the shadowing/ antishadowing effects in weak current 
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Figure 9: The nuclear shadowing and antishadowing effects at (Q 2 ) 
experimental data are taken from Refs. [4T| I48j . 



1 GeV 2 . The 



interactions. 



Nuclear effects on extraction of sin 2 6w 



The observables measured in neutrino DIS experiments are the ratios of neutral cur- 
rent (NC) to charged current (CC) current events; these are related via Monte Carlo 
simulations to sin 2 6w- In order to examine the possible impact of nuclear shadowing 
and antishadowing corrections on the extraction of sin 2 ^, one is usually interested 
in the following ratios 



R A 



R u a 



a{yp + A -> [T + xy 
(x(z7 M + A -> /2+ + X) 



(38) 
(39) 



of NC to CC neutrino (anti-neutrino) cross sections for a nuclear target A. As is well 
known, if nuclear effects are neglected for an isoscalar target, one can extract the 
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Figure 10: Our prediction for the ratio of F2 structure functions 
F 2 A [neutrino] /( 18/5) FaAfmuon], measured in Ref. jlO], at Q 2 = 20 GeV 2 . 



weak-mixing angle by using the Llewellyn-Smith relation [51] : 

„ u[u] _ gMgg] + N -> Vp] +X) ./I , , 5 ^4. n ,J-1]^ f4f)1 

cr(^[^J + N -> /i [/i+J+X) V2 9 / 

written in terms of NC and CC (anti-)neutrino-nucleon cross sections. Here, 

W cos 2 W M1' <r(i/ M + iV -> //- + X) 2' 1 j 

However, actual targets such as the iron target of the NuTeV experiment are not 
always isoscalar, having a significant neutron excess. In addition, as we have stressed 
here, nuclear effects due to multi-scattering could be very important. These nuclear 
effects should also modify the CC and NC structure functions, and therefore a detailed 
study of these effects on the extraction of the weak-mixing angle is essential. In order 
to reduce the uncertainties related to sea quarks, Paschos and Wolfenstein [H2] showed 
that one can extract sin 2 6w from the relationship 

a(u ll + N^Up + X)-<T(v ll + N^V li + X) 2 /l . 2 
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Inspired by the above relation, we will examine nuclear effects on sin 2 9w by the 
following observable for the scattering off a nuclear target A, 

o- °K + A^u^ + X)-a{V^ + A^V^ + X) 
A a^ + A^ ix- + X)-a(T7„ + A^ + 




Figure 11: Ratios F^/F^° (solid curves) and F^/F^° (dashed curves) for various 
current exchange interactions, at Q 2 = 1 GeV 2 . 

In the previous section, we have shown in Fig. El the nuclear effect on the electro- 
magnetic structure functions. Here we can also look the nuclear effect on the cross 
sections in CC and NC neutrino-nucleus DIS. In Fig. ^2 we show ratios F^/FJ?° 
(solid curves) and F^/Ff (dashed curves) for various current exchange interactions. 
The fact that the F A /F^° ratio for the VT~-current becomes negative and divergent 
for small x comes from the behavior of F^° , which in our model vanishes for x ~ 0.01. 
In addition, we are interested in the following ratios 

, , _ da{y„ + A-*v ll + X)/dx 
z[ ) da^ + N ^v lM + X)/dx 1 1 } 
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doiy^ + A — > + X) / dx 



(45) 




da(V^ + N V ll + X)/dx } 
daiv^ + A —>■ fi~ + X) / dx 



(46) 



daiv^ + N — > /i + X)/dx' 




(47) 



if 



'A/N 



R A /R 



N 



(48) 



which are later used to extract the nuclear effect on the weak-mixing angle. In Fig. 1121 
we show the ratios of Eqs. (|4*4"|) - (|4Tj) . From Figs. IHHT2*| one finds that the nuclear 
effect for charged and neutral currents is substantially different from that for the 
electromagnetic nuclear structure functions. There is a strong antishadowing effect 
in R v z and -R^-, but there is a small one in R u z and if^+. Moreover, for neutrinos 
the NC and CC shadowing/ antishadowing effects are the same, but for antineutrinos 
they are substantially different. As a result, in the neutrino case the Llewellyn-Smith 
relation can be used in order to extract the weak mixing angle, but it cannot be used 
for this purpose in the case of antineutrino deep inelastic scattering in nuclei. 

If the nuclear target had zero isospin, and there was no contribution from s or s 
quarks, then there would be no nuclear correction to the Llewellyn Smith relation. 
However, when one includes the strange quark currents, the situation is different. The 
neutral current interactions of the antineutrino on the s are shadowed rather than 
antishadowed (see Fig. ITTj) in the region x ~ 0.1, which reduces the total antishadow- 
ing effect. The charged current interactions of the antineutrino still experience strong 
antishadowing from the d. In contrast, neutrino interactions are relatively insensitive 
to the strange quark current since they are dominated by interactions on the valence 
quarks, not the anti-quarks. Even when one includes the strange quarks, antineutrino 
NC and CC interactions experience more antishadowing than neutrinos (see Fig. 112)) . 

An alternative way of assessing the nuclear corrections is through a modified 
Paschos-Wolfenstein ratio R^(x)/R]^(x), in which instead of the total cross section 
we consider the corresponding differential cross sections, as shown in Fig. ED In 
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Figure 12: The nuclear effect on the cross sections of CC and NC neutrino- nucleus 
DIS, at Q 2 = 1 GeV 2 . The dotted and dashed curves almost overlap. 
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Figure 13: The nuclear effect on the Paschos-Wolfenstein ratio of differential cross 
sections Rj(x)/R^(x), at Q 2 = 20 GeV 2 . 
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this case we have taken Q 2 = 20 GeV 2 , which is the average value of the NuTeV 
experiment. 

In our numerical analysis we studied the influence of nuclear effects on the extrac- 
tion of sin 2 9w from the observable R A , taking into account some kinematical cut-offs 
specific to the NuTeV experiment. 

The differential cross sections for CC and NC (anti-)neutrino-nucleus deep inelas- 
tic scattering are given by [S3] 

Gl m N E u>v {xy 2 F? ±{A \x,Q 2 ) + (49) 



dxdy 7T 



+ ii-»-^)Fr w (x,o»)±(»-^)«Fr w (*,^ 



for the CC reaction, and 



^pf A) = ^m N E u Jxy 2 F 1 Z{A \x,Q 2 )+ (50) 
dxdy 7r L 

+ (l - V - pg*) ff*><«, Q 2 ) ±(y-£\ ,f™<«,<«} , 



for the NC reaction. 

In the event selection, the NuTeV Collaboration applied the cut off 

20 GeV < E cal < 180 GeV, (51) 

for a visible energy deposit to the calorimeter E ca \. The lower limit ensures full 
efficiency of the trigger, allows for an accurate vertex determination and reduces 
cosmic ray background. 

Therefore we will calculate the observables R 1 ^^ and R A imposing the same cut 
off on the energy Eh of the final hadronic state X assuming Eh = E ca \. Since Eh ~ f, 
we can write the kinematical variables averaged over the (anti-)neutrino flux as 

Q2 < - Ecal < 

2M N E cal ~ ' V ~ (E u(v) ) ~ ' (52) 

For the average energies of the neutrino and antineutrino beams we take the values 
(E v ) = 120 GeV and (Ep) = 112 GeV, as in the NuTeV experiment 
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We will assume a modified version of the Paschos-Wolfenstein relation: 

1 



-R^sin 2 9 W ) 



1 



sin 2 # 



(53) 



P~ 



— sin 2 9 W 



where, p 2 = p%(l + e) with a modified factor (1 + e) due to strange quark, isospin 
breaking, threshold corrections for heavy quarks production, and so on. We further 
assume that the Paschos-Wolfenstein relation can be applied to the scattering on a 
nuclear target A, 



RJsin 2 9 W ) 



p 2 ( - - (sin 2 9 W + A sin 2 9 W ] 



(54) 



with a correction Asm 2 9w to the weak- mixing angle. In Fig. lT4T a). we show the 



sin 2 9w dependence in the ratio R A j N - We estimate A sin 2 6 W in the following way. 

2 , 



First, we use the cross sections to calculate the Paschos-Wolfenstein ratios R A (sm 9yy) 
and R N (sm 2 9w) with various values of sin 2 6^. Second, we extract p 2 by means of 
Eq. (J53)) . In principle, p 2 should be different for various values of sin 2 ^. We find a 
weak dependence of p 2 on sin 2 9w and p 2 ~ 1.04. Finally, we use the obtained p 2 to 
extract the shadowing/antishadowing effect on the weak-mixing angle A sin 2 9w from 
Eq. (}54j) . The results are given in Fig. ITiTb). 
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Figure 14: (a) The sin 2 9w dependence in R A / N \ 
ing/antishadowing corrections to the sin 2 9w 
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(b) The nuclear shadow- 
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We have performed a numerical calculation at the average (Q 2 ) = 20 GeV 2 of 
the NuTeV experiment and have found that the modification to the weak-mixing 
angle is approximately 5 sin 2 9w = 0.001. The value of sin 2 9w determined from the 
NuTeV experiment, without including nuclear shadowing/ ant ishadowing due to mul- 
tiple scattering, is in absolute value 0.005 larger than the best value obtain from other 
experiments. The model used here to compute nuclear shadowing/antishadowing ef- 
fect would reduce the discrepancy between the neutrino and electromagnetic mea- 
surements of sin 2 9]y by about 20%. Together with the charge symmetry violation 
contributions to the neutrino reactions [S3] , about half of the difference between the 
standard model and the NuTeV result can be accounted for. We also note that the 
antishadowing effects we predict are most important in the antineutrino data, which 
is less sensitive to sin 2 6 W . 

4 Conclusions 

We have investigated nuclear shadowing and antishadowing effects arising from the 
multiple scattering of quarks and anti-quarks in the nucleus. The effective quark- 
nucleon scattering amplitude includes Pomeron and Odderon contributions from multi- 
gluon exchange as well as Reggeon quark-exchange contributions. The model is con- 
strained by measurements of the nuclear structure functions in deep inelastic electron 
and muon scattering as well as the Regge behavior of the non-singlet structure func- 
tions. We have also noted the possibility of obtaining an antishadowing contribution 
from one-gluon exchange x Reggeon exchange, assuming the existence of hidden- 
color components in the nuclear wavefunction. We have shown that the coherence of 
these multiscattering nuclear processes leads to shadowing and antishadowing of the 
electromagnetic nuclear structure functions in agreement with measurements. The 
momentum sum rule is not satisfied in a nuclear target by balancing the shadowing 
and antishadowing of the leading-twist nuclear quark distributions; however the mo- 
mentum sum rule can still be satisfied if there is a compensating change in the nuclear 
gluon distribution. 

Our analysis leads to substantially different nuclear antishadowing for charged 
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and neutral current reactions; in fact, the neutrino and antineutrino DIS cross sec- 
tions are each modified in different ways due to the various allowed Regge exchanges. 
The non-universality of nuclear effects will modify the extraction of the weak-mixing 
angle sin 2 9w, particularly because of the strong nuclear effects for the F 3 structure 
function. The shadowing and antishadowing of the strange quark structure func- 
tion in the nucleus can also be considerably different than that of the light quarks. 
We thus find that part of the anomalous NuTeV result for sin 2 9w could be due to 
the nonuniversality of nuclear antishadowing for charged and neutral currents. Our 
picture also implies non-universality for the nuclear modifications of spin-dependent 
structure functions. 

We have found in our analysis that the antishadowing of nuclear structure func- 
tions depends in detail on quark flavor. Careful measurements of the nuclear de- 
pendence of charged, neutral, and electromagnetic DIS processes are thus needed 
to establish the distinctive phenomenology of shadowing and antishadowing and to 
make the NuTeV results definitive. It is also important to map out the shadowing 
and antishadowing of each quark component of the nuclear structure functions to 
illuminate the underlying QCD mechanisms. Such studies can be carried out in semi- 
inclusive deep inelastic scattering for the electromagnetic current at Hermes and at 
Jefferson Laboratory by tagging the flavor of the current quark or by using pion and 
kaon-induced Drell-Yan reactions. A new determination of sin 2 $w is also expected 
from the neutrino scattering experiment NOMAD at CERN [H7j . A systematic pro- 
gram of measurements of the nuclear effects in charged and neutral current reactions 
could also be carried out in high energy electron-nucleus colliders such as HERA and 
eRHIC, or by using high intensity neutrino beams 58,. 

Appendix: A non-abelian source for antishadowing 

We can identify a further antishadowing contribution specific to the non-abelian 
theory. Consider once again Fig. [2] for 7* A — > 7* A, but replace the two exchanged 
gluons with just a single gluon (see Fig. ITK|) . For simplicity we display the case 
of a deuteron target. In Fig. 4(a) the exchanged gluon attaches to the struck u 
quark valence constituent of the proton at the top of the diagram changing its color. 
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Figure 15: Representation of leading-twist QCD "hidden color" contributions to the 
nuclear structure function from the absorptive part of the forward virtual Compton 
amplitude Im T{<y*A — > 7* A), (a) Illustration of a two-step contribution in the q + = 
0) Q± = Q 2 parton model frame — deep inelastic lepton scattering on a valence quark 
of a target proton followed by the final-state single-gluon interaction of the outgoing 
quark on a target neutron. The proton and neutron are both color excited to color- 
octet states. The amplitude requires the presence of hidden-color components in 
the nuclear wavefunction. (b) Illustration of the physics of the two-step process 
shown in (a), but in the laboratory frame where q + > 0. The uu fluctuation of the 
virtual photon first scatters via a single gluon exchange on a neutron; this is then 
followed by the annihilation of the u quark on a proton. The proton and neutron 
are both color excited to color-octet states. The amplitude requires the presence of 
hidden-color components in the nuclear wavefunction. The two-step amplitude can 
interfere constructively with the Regge-behaved single-step annihilation amplitude on 
the proton alone, thus producing antishadowing. Similar processes occur in the case 
of the weak currents. 
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This also changes the scattered proton p' to a color octet. The exchanged gluon also 
transforms the spectator neutron into a color octet. Thus if the deuteron wavefunction 
contains hidden color |8c8c) components, this process interferes with the one step 
diagram with no final state interactions. 

The deuteron certainly has hidden-color components — one only has to exchange a 
gluon between the nucleons in the deuteron LFWF The large magnitude of the 
deuteron form factor also demands hidden color components ([HE])- The calculation 
of the one-gluon exchange effects is very similar to our Odderon analysis. The one- 
gluon exchange amplitude behaves as s 1 and a nearly real phase. Like the Odderon, 
it has C = — and couples with opposite sign to the q and q. A complication is how to 
understand the Reggeon exchange amplitude on the proton since the u and u in the t 
channel now are in a color octet configuration. Nevertheless, the Odderon calculation 
serves as a model for the one-gluon exchange contribution as well and its effect on 
antishadowing. 
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